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Abstract 

Geoengineering with stratospheric sulfate aerosols has been proposed as a means of temporarily 
cooling the planet, alleviating some of the side effects of anthropogenic CO 2 emissions. However, one 
of the kn own side effects of stratospheric injections of sulfate aerosols is a decrease in stratospheric 
ozone. Here we show results from two general circulation models and two coupled chemistry climate 
models that have simulated stratospheric sulfate aerosol geoengineering as part of the Geoengineering 
Model Intercomparison Project (GeoMIP). Changes in photolysis rates and upwelling of ozone -poor air 


in the tropics reduce stratospheric ozone, suppression of the NO x cycle increases stratospheric ozone, 
and an increase in available surfaces for heterogeneous chemistry modulates reductions in ozone. On 
average, the models show a factor 20-40 increase of the sulfate aerosol surface area density (SAD) at 
50 hPa in the tropics with respect to unperturbed background conditions and a factor 3-10 increase at 
mid-high latitudes. The net effect for a tropical injection rate of 5 Tg SO 2 per year is a decrease in 
globally averaged ozone by 1.1 -2.1 DU in the years 2040-2050 for three models which include 
heterogeneous chemistry on the sulfate aerosol surfaces. GISS-E2-R, a fully coupled general 
circulation model, performed simulations with no heterogeneous chemistry and a smaller aerosol size; 
it showed a decrease in ozone by 9.7 DU. After the year 2050, suppression of the NO x cycle becomes 
more important than destruction of ozone by C10 x , causing an increase in total stratospheric ozone. 
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Contribution of ozone changes in this experiment to radiative forcing is 0.23 W m' in GISS-E2-R and 
less than 0.1 W m" in the other three models. Polar ozone depletion, due to enhanced formation of both 
sulfate aerosol SAD and polar stratospheric clouds, results in an average 5% increase in calculated 
surface UV-B. 


1. Introduction 

Geoengineering has been proposed as a means of deliberately, temporarily altering the climate to 
alleviate some of the consequences of anthropogenic greenhouse gas emissions [e.g., Shepherd et al . , 
2009]. One proposed method involves mimicking a volcanic eruption by injecting sulfate aerosol 
precursors (e.g., SO 2 ) into the stratosphere [ Budyko , 1974; Crutzen, 2006]. Large volcanic eruptions 
are capable of injecting considerable amounts of particles and sulfur gases above the tropopause, 
causing increases in stratospheric aerosol optical depth, sometimes over one order of magnitude. These 
volcanic events, the largest most recent ones being Agung, El Chichon, and Pinatubo, in 1963, 1982, 
and 1991, respectively, can cause substantial surface cooling. The 1991 eruption of Pinatubo, for 
instance, resulted in reductions in globally averaged surface air temperature by approximately 0.5°C 
[Soden et al., 2002]. Although volcanic aerosols are episodic, whereas many proposed methods of 
geoengineering involve sustained injections, they provide important clues as to some of the expected 
climate response to geo engineering. 

Chemical-radiative-dynamical perturbations of the stratosphere following the 1991 eruption of Mt. 
Pinatubo have been widely studied by the scientific community and reported in the literature. 


Nucleation of H 2 SO 4 vapour derived from the initial SO 2 plume [Bluth et al, 1992; Read et al., 1993] 
had produced an optically thick cloud of sulfate aerosols [ McCormick and Veiga, 1992; Lambert et al, 
1993; Long and Stowe, 1994], which are highly reflective in the visible and UV. This effect is 
particularly important in the tropics, where the aerosols (depending on the QBO phase) may remain 
confined for several months after the eruption with relatively high values of optical thickness, as it was 
the case after Pinatubo [ Trepte and Hitchman, 1992]. The radiative feedback of volcanic aerosols is 
due to additional diabatic heating introduced by: (a) direct interaction with solar and planetary radiation 
and indirect ozone absorption of the aerosol diffused solar radiation \Pitari, 1993; McCormick et al., 
1995]; and (b) changing ozone absorption of incoming solar and outgoing planetary radiation, as a 
consequence of ozone losses produced by heterogeneous chemistry on volcanic aerosols, perturbed 
photolysis rates and stratospheric large scale transport [Prather, 1992; Kinne et al., 1992; Schoeberl et 
al, 1993]. Perturbations of stratospheric tracer species transport may result from dynamical changes 
due to both local stratospheric heating and climate changes associated to the increasing scattering of 
incoming solar radiation by the volcanic particles [Hansen et al, 1992; Lads et al., 1992; Stenchikov et 
al., 1998; Kirchner et al, 1999]. The stratospheric mean meridional circulation is affected by local 
aerosol radiative heating (mostly located in the tropical lower stratosphere). The planetary wave 
propagation in the mid- to high-latitude lower stratosphere is altered as a consequence of changing 
atmospheric stability due to the climate perturbation. The 1991-1992 time series of the calculated 
globally averaged net radiative flux change at the tropopause due to aerosols from the Pinatubo 
eruption [Hansen et al, 1992] shows that the largest forcing occurs during January 1992 (-4.5 W/nr, 
with x = 0.15); the magnitude of radiative flux changes decreases afterwards with an e-folding time of 
about one year (-1.2 W/m during June 1993). Measurements taken a few months after the Pinatubo 
eruption revealed a 2-3 K wanning in the tropical lower stratosphere [Labitzke and McCormick, 1992; 
Young et al., 1994] and a 20-25% ozone depletion, with a decrease of about 20 DU of the ozone 
column in the 16-28 km layer in the tropics during October-November 1991 [Grant et al., 1992; 
Schoeberl et al., 1993]. 

In summary, the major effects on stratospheric ozone of large volcanic perturbations of stratospheric 
aerosols are obtained through: (a) photolysis (of O2 and O3, in particular); (b) heterogeneous chemistry 
on the surface of sulfuric acid aerosols; (c) homogenous chemistry, via temperature changes; (d) 
heterogeneous chemistry on the surface of polar stratospheric clouds (PSC), via temperature changes 
and increasing population of sulfate aerosols condensation nuclei; and (e) perturbations of the 


stratospheric large scale circulation. Photochemical processes (a)-(c) are mostly relevant for the tropics 
and mid-latitudes, whereas effect (d) could be important for polar ozone depletion. Although the 
stratospheric effects due to volcanic eruptions are well studied, the expected effects due to 
geoengineering have received less attention. Tilmes et al. [2008] reported that geoengineering with 
stratospheric sulfate aerosols could enhance stratospheric ozone destruction and delay the recovery of 
the Antarctic ozone hole by 30-70 years. Further, the acceleration of ozone loss cycles over the polar 
regions due to the increase of aerosols may significantly decrease the ozone column and increase 
erythemal UV by up to 5% in mid and high latitudes and 10% over Antarctica, by mid-century halogen 
conditions [ Tilmes et al., 2012]. 

The Geoengineering Model Intercomparison Project (GeoMIP) was initiated to determine the robust 
features of climate model response to four core geoengineering experiments involving uniform solar 
reduction and stratospheric injections of sulfate aerosol precursors [Kravitz et al, 2011a]. In this study, 
we investigate the stratospheric ozone response to experiments G3 and G4, with focus on 
photochemical perturbations. Both of these experiments have a background anthropogenic forcing 
profile corresponding to RCP4.5 [ Taylor et al., 2012]. Beginning in 2020, G3 involves transient 
stratospheric injections of SCfi to maintain top of atmosphere net radiation at 2020 levels. G4 involves 
a constant stratospheric injection rate of 5 Tg SO 2 per year. These simulations are performed for 50 
years, after which geoengineering is immediately ceased; the simulations are then run for an additional 
20 years to determine the climate response of the termination effect [e.g., Jones et al., 2013]. Because 
all participating models have performed these same experiments, we can determine the climate model 
response to these scenarios independent of inter-model differences. In particular, we can determine 
robust ozone responses in climate models to a sustained layer of stratospheric aerosols 
The paper is organized in three sections. Section 2 includes a description of participating models. 
Section 3 discusses models results regarding stratospheric aerosols (background and GeoMIP 
perturbations). Section 4 analyzes changes in stratospheric ozone in terms of aerosol induced 
perturbations in chemical production and loss, vertical profiles, total column, tropopause radiative 
forcing and surface UV. The main conclusions are summarized in section 5. 


2. Description of models 

Participating models in this study are the following: ULAQ-CCM (University of L’Aquila Chemistry 
Climate Model), GISS-E2-R (Goddard Institute for Space Science Model E2), MIROC-ESM-CHEM 


(Model for Interdisciplinary Research on Climate), GEOSCCM (Goddard Earth Observing System 
Chemistry Climate Model). A summary of model runs is presented in Table 1. 

The main differences in the treatment of stratospheric aerosols and ozone are summarized here, (a) 
Aerosol formation starting from oxidation of SO 2 is included in ULAQ-CCM, GISS-E2-R, 
GEOSCCM; explicit prediction of the aerosol size distribution from microphysical processes is made 
in ULAQ-CCM; a fixed lognormal size distribution for the dry sulfate aerosol, hydrated depending on 
the relative humidity, is present in GEOSCCM; MIROC-ESM-CHEM does not calculate aerosol 
fonnation in the stratosphere, instead aerosol optical depth (AOD) is specified as a function of latitude, 
height, month with a fixed lognormal size distribution used in radiative calculations (0.237 |un 
effective radius); a fixed gamma size distribution for the dry sulfate aerosol, hydrated depending on 
relative humidity, is used in GISS-E2-R. Aerosol heating rates in solar near-infrared (NIR) and 
planetary longwave radiation are included in all models, (b) Ozone photochemistry with sensitivity of 
photolysis to aerosol scattering is in ULAQ-CCM, GISS-E2-R, MIROC-ESM-CHEM; a look-up table 
for photo-dissociation rates is used in GEOSCCM. (c) Heterogeneous chemical reactions on the surface 
of sulfate aerosols, with feedback on NO x chemistry, is included in ULAQ-CCM, MIROC-ESM- 
CHEM, GEOSCCM; it is not included in GISS-E2-R and in specific sensitivity runs of ULAQ-CCM . 
The acronym ‘nhc’ used throughout the paper stands for ‘no heterogeneous chemistry on stratospheric 
sulfate aerosol’. 


2.1 ULAQ-CCM 

The University of L’Aquila model is a low-resolution global scale climate-chemistry coupled model 
(ULAQ-CCM) extending from the surface to the mesosphere (0.04 hPa); model features are described 
by Pitari et al. [2002], Eyring et al. [2006], and Morgenstern et al. [2010]. After participation in 
CCMVal-2, some important updates have been made to the model: (a) increase of horizontal and 
vertical resolution, now T21 with 126 log-pressure levels (approximate pressure altitude increment of 
568 m); (b) inclusion of a parameterization for the formation of upper tropospheric cirrus cloud ice 
particles [Karcher and Lohmann, 2002]; (c) update of species cross sections using recommendations by 
Sander et al. [2011] and Schumann-Runge bands treatment following the parameterization of 
Minschwaner et al. [1993] based on fixed-temperature opacity distribution function formulation; (d) 
upgrade of the radiative transfer code for calculations of photolysis, solar heating rates and top-of- 


atmosphere radiative forcing. The oceanic surface temperature is assimilated from the Hadley Centre 
for Climatic Prediction and Research [Rayner et al., 2003]; the parameterization of periodic natural 
forcings (solar cycle, quasi-biennial oscillation) is included on-line. The chemistry module is organized 
with all medium and short-lived species grouped in the families O x , NO y , NO x , CHO x , Cl y , Br y , SO x , 
and the module also includes aerosols. In total there are 40 transported species, 26 species at 
photochemical equilibrium and 57 size categories for aerosols. The model includes the major 
components of stratospheric and tropospheric aerosols (sulfate, carbonaceous, soil dust, sea salt, PSCs). 
The size distributions of sulfate (both tropospheric and stratospheric) and PSC particles are calculated 
using an interactive and mass conserving microphysical code for aerosol formation and growth, 
including a gas-particle conversion scheme. New sulfuric acid particles are formed via homogeneous 
nucleation, followed by coagulation and condensation growth; stratospheric aerosols are lost via 
evaporation in the upper stratosphere and downward transport and sedimentation into the troposphere. 
For the GeoMIP simulations, SO 2 was injected throughout the altitude range of 18-24 km at 0° 
longitude on the equator. Lower stratospheric denitrification and dehydration are calculated using the 
predicted size distribution of PSC particles; heterogeneous reactions on sulfate and PSC aerosols are 
included in the chemistry module [ Pitari et al., 2002], with the hydrolysis of N 2 0 5 and BrONCL being 
the two most important reactions on the surface of stratospheric sulfuric acid aerosols . The updated 
radiative transfer module, operating on-line in the ULAQ-CCM, is a two-stream delta-Eddington 
approximation model [ Toon et al, 1989], used for chemical species photolysis rate calculation in ultra- 
violet (UV)-visible (VIS) wavelengths and for solar heating rates and radiative forcing in UV-VIS- 
near-infrared (NIR) bands. Top-of-atmosphere solar fluxes are taken from SUSIM-SL2 and 
LOWTRAN7 and are integrated on the wavelength bins used in the model, of which there are 150 in 
the UV and visible range and 100 in the NIR range, covering the solar spectrum from Lyman-a up to 7 
pm. Sun-earth distance is calculated daily as a function of orbit eccentricity and the solar cycle is 
included. Sphericity is treated by means of Chapman functions [Dahl back and Stamnes, 1993]. 
Refraction is taken into account with an iterated ray-tracing technique in a simple exponential 
refraction model. Optical depths take into account Rayleigh scattering, radiation absorption from O3, 
0 2 , NO 2 , SO 2 , H 2 0, C0 2 and scattering/absorption from aerosol particles. Planetary radiation heating 
rates and top of atmosphere/tropopause forcings are calculated including absorption/emission by CO 2 
(15 pm band), O3 (9.6 pm band), H 2 0 and aerosols, using well tested absorptivity fonnulas through a 
correlated-k method. Aerosol extinction values are passed daily from the ULAQ-CCM aerosol module 
to the radiative transfer module, with appropriate wavelength-dependent values of extinction 


coefficient, asymmetry factor and single scattering albedo, given the calculated size distribution of the 
particles. Surface albedo is taken from MERRA 2D hourly averaged data. The ULAQ model 
calculations of photolysis rates and surface and top-of-atmosphere radiative fluxes have been validated 
in the framework of CCMVal and AEROCOM inter-comparison campaigns [ Chipperfield and Liang et 
al, 2013; Randles et al., 2013]. 


2.2 GISS-E2-R 

GISS-E2-R is a fully coupled atmosphere-ocean general circulation model, developed by the National 
Aeronautics and Space Administration’s (NASA) Goddard Institute for Space Studies (GISS), which is 
a contributor to the Coupled Model Intercomparison Project Phase 5 (CMIP5) [ Schmidt et al., 2006; 
Taylor et al., 2012]. The atmospheric model has a horizontal resolution of 2° latitude by 2.5° longitude 
and 40 vertical layers extending through the mesosphere (model top of 0.1 hPa). It is coupled to the 
Russell ocean model [ Russell et al., 1995] which has horizontal resolution of 1° latitude by 1.25° 
longitude and 32 vertical layers. The version discussed here was run with full stratospheric chemistry, 
involving 113 fully resolved reactions. Photolysis in the stratosphere is included and varies 
dynamically as a function of temperature, pressure, solar output, cloudiness, albedo, and ozone 
concentration. Ozone and methane are computed at the model’s resolution and are radiatively 
interactive with model tracers, i.e., they modify heating rates and are chemically active. Formation of 
sulfate aerosols is specified by a reaction rate in which an aerosol dry radius is specified, as described 
by Koch et al. [2006]. The aerosols then grow hygroscopically according to ambient relative humidity 
per the formulas of Tang [1996]. In this study, the aerosol dry radius is specified to be 0.35 pm, which 
is the same specification as was used by Robock et al. [2008] in their simulations of the 1991 eruption 
of Mount Pinatubo and results in an approximately Pinatubo-sized aerosol in simulations of volcanic 
eruptions, per GeoMIP specifications [Kravitz et al., 2011b]. SCF was injected throughout the altitude 
range of 16-25 km at 0° longitude on the equator. The aerosols are radiatively active tracers and are 
advected within the model via the general circulation. Stratospheric sulfate aerosols injected into the 
tropical lower stratosphere have an e-folding lifetime of approximately 12 months. 


2.3 MIROC-ESM-CHEM 


MIROC-ESM-CHEM is an Earth System Model (ESM) that has been developed based on a global 
climate model MIROC (Model for Interdisciplinary Research on Climate) consisting of coupled 
atmosphere, ocean, sea ice, river, and land surface models. The details of MIROC-ESM-CHEM model 
and the base (RCP4.5) simulation settings have been documented in Watanabe et al. [2011a] and 
Watanabe et al. [2011b], respectively. The atmosphere model of MIROC-ESM-CHEM has a T42 
horizontal resolution (approximately 300 km grid spacing), and contains 80 vertical layers from the 
surface to a height of about 85 km. The model predicts effects of the solar cycle, QBO, and ENSO 
online. The model considers the direct and indirect effects of tropospheric aerosols, while treatments of 
the stratospheric aerosols are separated from the tropospheric one as is outlined below. 

In the base simulation, a zonally averaged stratospheric aerosol AOD is provided to the model's 
radiation module as a function of latitude, altitude and month following Sato et al. [1993], which is 
exponentially reduced after 1998 with one year relaxation time toward a background value of 10' 4 . In 
G4, we used similar latitude- altitude- and monthly data of AOD but provided by GeoMIP for models 
which do not simulate formation of the stratospheric sulfate aerosols from SO 2 gas \Kravitz et al., 
2011b]. hi radiation calculations, a constant effective radius of 0.237 pm is used in both the base and 
G4 runs. More details of the radiation calculations are discussed by Watanabe et al. [2011b]. 

The atmospheric chemistry module of MIROC-ESM-CHEM predicts the system of tropospheric 
chemistry (O x -HO x -NO x -CO-CH 4 - VOCs: Volatile Organic Compounds) as well as major chlorine and 
bromine compounds (Cl y and Br y ) that are important for the simulation of stratospheric ozone 
[Watanabe et al., 2011a]. Parameterizations for liquid and solid particles in the stratosphere are 
included to calculate heterogeneous reactions on liquid sulfate aerosols (including the hydrolysis of 
N 2 O 5 and BrONCh) and PSCs [e.g., Carslaw et al, 1995; Hanson and Mauersberger, 1988]. SADs of 
PSCs are diagnosed as a function of the predicted temperature and mixing ratio of H 2 O and HNO 3 . On 
the other hand, SAD of liquid sulfate aerosols ( S) in the stratosphere including those background and 
volcanic aerosols are given as an external forcing in this study. Namely, S (in pm"cm' ) is 
approximately diagnosed in the heterogeneous chemistry package based on the AOD data: 
Ar 1 

S = 4c , where c is a constant for unit conversion ( 1 0 10 ), Ax is x in a certain model layer at a 

A Z £ 

certain location, Az is a thic kn ess of the model layer (m), and s is a prescribed constant representing 
mean efficiency of scattering due to the sulfate aerosols (assumed to be 2). This relationship is derived 
through a combination of several basic and approximated equations, and gives S distribution which 


resembles to observations [e.g., Thomason et al, 1997b], via solving out n and o from the following: 
S = Art ■ r 2 n , At = noAz , and cr = sn ■ r 2 , where n, r and o are the number density, effective radius, 
and a mean scattering cross section of sulfate aerosols in the model layer. 


2.4 GEOSCCM 

The Goddard Earth Observing System Chemistry Climate Model (GEOSCCM) used for the GeoMIP 
simulations couples the Goddard Earth Observing System, version 5 (GEOS-5) [ Rienecker et al., 2008] 
general circulation model to the Georgia Institute of Technology-Goddard Global Ozone Chemistry 
Aerosol Radiation and Transport (GOCART) module [Chin et al, 2000, Colarco et al, 2010] and a 
stratospheric chemistry module [Pawson et al., 2010]. 

GEOS-5 uses a finite volume dynamical core [Lin, 2004] combined with a physics package that 
describes moist processes, radiation, turbulent mixing and surface processes. Convection is 
parameterized using the relaxed Arakawa-Schubert (RAS) scheme [Moorthi and Suarez, 1992], and is 
combined with a prognostic cloud scheme. The boundary layer turbulent mixing is parameterized with 
the schemes by Louis et al. [1982] and Lock et al. [2000] for stable and unstable situations, 
respectively. The land surface model is composed of a catchment-based hydrological model [Koster et 
al., 2000] and a multilayer snow model [Stieglitz et al., 2001]. The radiative transfer model consists of 
a solar radiation model [Chou and Suarez, 1999] and a thennal radiation model [Chou et al., 2001]. 
The solar radiation model includes absorption due to water vapor, O3, O2, CO 2 , clouds and aerosols. 
The thermal radiation model includes absorption by water vapor, CO 2 , O3 and most of the minor trace 
gases, as well as clouds and aerosols. GOCART includes a parameterization of the chemical production 
of SO 4 aerosol from oxidation of dimethyl sulfide (DMS) by OH during day and NO3 during night, and 
from oxidation of SO 2 by OH in the gas phase and by H 2 O 2 in the aqueous phase. 

The surface area density of stratospheric aerosols is used in the heterogeneous chemistry module. 
GEOSCCM calculates the aerosol SAD from the concentration of SO 4 assuming that the dry 
stratospheric sulfate aerosol particles are lognonnally distributed with modal radius equal to 0.35 pm 
and standard deviation 1.6. These values have been chosen within the observed range [e.g. Bingen et 
al., 2004] and result in good agreement with observations of the period after the eruption of Mount 


Pinatubo [Aquila et al., 2012]. The dry stratospheric aerosol is hydrated depending on the ambient 
relative humidity following Petters and Kreidenweis [2007]. 

A similar version of GEOSCCM has been evaluated by Aquila et al. [2012] and Aquila et al. [2013]. 
Respect to Aquila et al. [2013], the version of GEOSCCM used in this work includes the coupling 
between aerosol and heterogeneous chemistry through the aerosol SAD, different parameters for the 
aerosol size distribution to reach a better agreement with observations, and an internal mechanism for 
the generation of the QBO [ Molod et al., 2012]. GEOSCCM is used in this work to perfonn only G4 
simulations and the relative control simulations. Here, 5 Tg per year of SO 2 is continuously injected at 
the equator at 0° longitude. The injection is uniform between 16 km and 25 km altitude. The 
transformation of SO 2 into sulfate aerosol is calculated by GOCART using climatological oxidant 
fields [Chin et al, 2000]. 

GEOSCCM does not directly simulates the stratospheric background aerosol. In the RCP4.5 control 
simulations, GEOSCCM prescribes aerosol SAD from the Stratospheric Aerosol Measurement II 
(SAM II) and Stratospheric Aerosol and Gas Experiment (SAGE) data ( Eyring et al., 2008) relative to 
the year 1979 - a period when the stratospheric aerosol layer was relatively unperturbed ( Thomason et 
al, 1997a)- for the purpose of calculating the heterogeneous chemistry. The radiative effect of the 
background aerosol, however, is not included in this simulation. In the G4 simulations the SAD 
calculated from the simulated sulfate aerosol mass is added to the prescribed background SAD. 

The model resolution used for these simulations is 2.0° latitude by 2.5° longitude and 72 vertical hybrid 
levels up to 0.01 hPa. The model is forced with sea surface temperatures and sea ice concentrations 
from RCP4.5 future simulations performed with the Community Earth System Model (CESM) [Gent et 
al, 2011]. 


3, Stratospheric aerosols 

Explosive volcanic perturbations and potential stratospheric sulfate aerosol geoengineering may act as 
large stratospheric sources of SO 2 , producing transient or steady state major increases of the 
stratospheric sulfate aerosol number, mass and surface area densities, as well as solar radiation 
extinction and optical depth; this can result in significant changes of the particle size distribution 
respect to the stratospheric background aerosol [Deshler et al., 1992; Thomason et al., 1997b]. These 


large aerosol perturbations directly impact the radiation budget at the surface and at the top of 
atmosphere, as well as the hydrologic cycle and tropospheric and stratospheric ozone [e.g., Robock, 
2000; Robock et a/., 2008; Tilmes et al, 2008]. To understand the stratospheric ozone anomaly 
induced by a perturbation of the stratospheric sulfate aerosol burden, we must first establish the aerosol 
size distributions and surface area. For example, larger aerosols will have a greater infrared radiative 
effect with respect to small aerosol, causing increased stratospheric heating and hence greater radiative- 
chemical-dynamical changes. A change in surface area density will change the amount of surface 
available for heterogeneous chemistry. 

Background stratospheric sulfate aerosols are either formed locally by OH oxidation of SO 2 or 
transported upward from the troposphere through the tropical tropopause layer (TTL), although the 
latter process is largely limited by upper tropospheric formation of cirrus ice particles via 
homogeneous freezing on sulfate aerosols [ Hendricks et al., 2011]. Stratospheric SO 2 , in turn, 
originates either from convection of SO 2 from the boundary layer to the TTL, followed by vertical 
advection from the TTL, or from local photochemical production, following photolysis of OCS [Pitari 
et al, 2002; Weisenstein et al., 2006]. Sulfuric acid produced at the end of the SO 2 oxidation chain 
forms aerosol particles via homogeneous nucleation (and to a lesser extent via heterogeneous 
nucleation of carbonaceous particles transported from the troposphere). The particle size distribution is 
detennined by these and other microphysical processes, such as condensational growth, coagulation, 
gravitational settling and evaporation in the upper stratosphere. The main sink for stratospheric sulfate 
aerosols is the mid-high latitude downward transport towards the troposphere coupled to gravitational 
sedimentation, with an average lifetime of about one year [ Weisenstein et al, 2006]. 

Model predictions of AOD perturbations in the GeoMIP experiments are presented in Fig. 1. G4 
imposes the amount of SO 2 injected mass flux (5 Tg/year): inter-model differences arise from longer or 
shorter stratospheric aerosol lifetimes, which in turn are produced by inter-model differences of the 
tropical pipe isolation and net mixing with extratropics, as well as differences of the Brewer-Dobson 
circulation strength and its potential modulation by stratospheric aerosol heating and QBO. The aerosol 
particle size distribution also plays an important role, by affecting the gravitational sedimentation rate. 
On the other hand, G3 imposes the conservation of the top-of-atmosphere radiative forcing (TOARF) 
during the GeoMIP time span (2020-2070). Assuming that the net change of the TOARF in these 50 
years is of the order of 2 W/m [ Meinshausen et al., 2011], the imposed stratospheric mass flux of SO 2 
in G3 may have large variations among the models. This variation depends mostly on the treatment of 


the aerosol effects on other atmospheric processes (e.g., O 3 , tropospheric cloudiness, albedo), since the 
direct aerosol equivalent-albedo effect is only slightly affected by the model-predicted aerosol size 
[Lacis et al., 1992]. 

In the present study, G3 with interactive ozone is simulated by two models (ULAQ-CCM and GISS- 
E2-R), both of which include explicit aerosol feedbacks with ozone photochemistry but not aerosol- 
cloud interactions. Model predicted aerosol optical depth for G3 is approximately a factor of 5 larger 
in ULAQ-CCM than GISS-E2-R. This is in part due to a low responsiveness of GISS-E2-R to CO 2 
changes as compared to several other models participating in CMIP5 [e.g., Kravitz et al., 2013], and 
consequently to a lower TOARF to be compensated by the geoengineering aerosol. The two models 
agree quite well in the predicted globally-averaged AOD in G4, where the SCE mass flux is imposed, 
independent of its net effect on the radiative forcing. 

Panels (d, e) of Fig. 1 show that the AOD latitudinal dependence is rather different among models that 
predict aerosols from the SO 2 injection: GEOSCCM and ULAQ-CCM have similar aerosol effective 
radii (Table 2) and then comparable stratospheric aerosol loss due to sedimentation. This, coupled to a 
good tropical pipe isolation [e.g. Strahan et al, 2011; Chipperfield and Liang et al., 2013], produces an 
effective aerosol confinement in the tropics, with a clear equatorial AOD maximum and comparable 
globally-averaged AOD. On the other hand, the GISS-E2-R latitudinal distribution shows a much faster 
aerosol dispersion towards mid-high latitudes (to be associated to a more efficient sub-tropical 
horizontal air mixing in the lower stratosphere), thus producing a tropical minimum of the AOD. The 
globally-averaged value, however, is comparable and even higher than that of GEOSCCM and ULAQ- 
CCM, due to the slower sedimentation loss produced by the smaller particle size (Table 2). 

The growth of stratospheric aerosol particles, produced by both transient or sustained SO 2 injections, 
modifies the particle size distribution shape leading to significantly larger effective radii with respect to 
normal background conditions of the stratosphere. This is highlighted both in direct satellite 
measurements of aerosol extinction after a major volcanic eruption [ Thomason et al., 1997b] and in 
global modeling studies [ Weisenstein et al., 2006; Heckendorn et al, 2009]. Owing to their large size, 
these particles have a reduced albedo and a faster sedimentation rate, with a net reduction of the 
stratospheric lifetime: both effects end up limiting the cooling potential with respect to stratospheric 
background aerosols. In addition, the increasing particle size favors the absorption of infrared radiation 
both in the solar NIR and planetary spectra, leading to larger additional tropical heating rates. This 


local radiative perturbation may have a potential increasing impact on tropical upwelling and (by 
continuity) on the extra-tropical downwelling, acting again to reduce the lifetime of stratospheric 
aerosols. Table 2 shows that two of the three models that form aerosols from the stratospheric SO2 
injection actually predict an increase of the effective radius in G4 with respect to RCP4.5 (i.e. ULAQ- 
CCM and GISS-E2-R), and two models calculate a G4 perturbed effective radius close to SAGE-II 
observations after the Pinatubo eruption (ULAQ-CCM and GEOSCCM). The effective radius is the 
radiative-effective aerosol dimension, detennining the efficiency of radiation scattering and absorption. 
Inter-model differences in Table 2 need to be taken into account to evaluate the GeoMIP aerosol impact 
on radiative processes, such as photolysis and stratospheric heating rates. 

The SAGE-II measured extinction ratio between visible and NIR wavelength channels (i.e. 0.525 and 
1.02 pm, respectively) is a good measure of the aerosol size distribution shape and its modifications 
produced by changes in the local sulfate productions rate from OH oxidation of SO2 [ Weisenstein et al., 
2006]. As shown in Fig. 2, this ratio is close to unity when the effective radius is close to 0.6 pm and 
about 3 for r e ff = 0.2 pm. The ULAQ-CCM is the only model included in this paper with explicit 
aerosol microphysics and calculates the time behavior of the sulfate aerosol size distribution as a 
function of stratospheric SO2, OH, and large scale transport. Its predictions in terms of the extinction 
ratio (and hence of the aerosol size distribution shape) are shown in Fig. 2 to be in good agreement with 
SAGE-II observations, for both background and perturbed conditions (G4 and Pinatubo). Panel (b) of 
Fig. 2 compares directly the ULAQ-CCM calculations of the aerosol size distribution for background 
and G4 conditions, confirming the previous findings. 

Aerosol SAD changes are presented in Fig. 3-4: these are important for assessing changes in 
heterogeneous chemical reactions in the stratosphere, particularly NO x and other O3 chemical 
precursors (Cl x , Br x ). The models use different approaches for calculating the SAD. ULAQ-CCM 
calculates directly the geometric surface area starting from the predicted aerosol size distribution. 
GEOSCCM calculates the SAD from the sulfate aerosol mass concentration, assuming that the dry 
stratospheric sulfate aerosol is log-normally distributed with modal radius 0.35 pm and standard 
deviation 1.59. This dry size distribution is hydrated depending on the relative humidity following 
Petters and Kreidenweis [2007]. MIROC-ESM-CHEM approximately calculates the SAD from the 
prescribed AOD distribution, as mentioned in Section 2.3. Model results are also validated with 
SAGE-II data, for both volcanically quiet and perturbed conditions. Although some differences are 
present in the latitudinal distribution of the aerosol SAD, all three models are well consistent with the 


reported SAGE-II derived SAD values for background and post-Pinatubo conditions. It should be noted 
that the average sulfate mass loading above the tropopause is comparable in G4 steady-state conditions 
(4.83 ± 0.06 Tg-S) with respect to April 1992 - March 1993 post-Pinatubo conditions (4.7 ± 0.6 Tg-S); 
these mass burden estimates are calculated as an average from ULAQ-CCM and GEOSCCM time- 
dependent simulations. This makes possible and meaningful to use the indirect G4 SAD validation with 
SAGE-II derived values from April 1992 to March 1993. On average, the three models show a factor 
20-40 increase of the tropical aerosol SAD at 50 hPa with respect to RCP4.5 (with ULAQ-CCM on the 
upper limit), and a factor 3-10 increase at mid-high latitudes (with GEOSCCM on the lower limit). 
These large sulfate aerosol SAD increases greatly affect the stratospheric NO x -HO x -Cl x -Br x chemistry 
[Fahey et al., 1993] by converting more NO x into HNO 3 and by increasing the amount of HO x and 
reactive Cl-Br. The balance of these catalytic cycle changes may enhance or limit the stratospheric 
ozone destruction depending on latitude, altitude and on the (time-dependent) total stratospheric 
amount of inorganic Cl y and Br y . 


4 . Stratospheric ozone 

The increase in stratospheric sulfate concentrations affects ozone via the enhancement of 
heterogeneous chemistry on the larger aerosol surface area density [e.g. Brasseur and Granier, 1992], 
changes in photolysis rates due to the increased aerosol optical thickness [ Kinne et al., 1992], and 
modifications of the atmospheric dynamics due to aerosol-ozone perturbations of the stratospheric 
heating rates [ Pitari and Rizi, 1993; Aquila et al., 2013], which are coupled to altered amplitude and 
propagation of planetary waves [Pitari and Mancini, 2002]. The temperature perturbations in the 
middle stratosphere may also affect the homogeneous chemical reaction rates with pronounced 
temperature dependence, such as O+O3 and NO+O3. ULAQ-CCM and MIROC-ESM-CHEM include 
all of these effects. GEOSCCM includes the ozone perturbations due to heterogeneous and 
homogeneous chemistry and to atmospheric dynamics, but not changes in photolysis rates. GISS-E2-R 
includes the ozone perturbations due to changes in photolysis rates, atmospheric dynamics and 
homogeneous chemistry, but not those due to heterogeneous chemistry on the sulfate aerosol surface. 
However, GISS-E2-R does include treatments of heterogeneous chemistry in polar stratospheric 
clouds, which is particularly relevant for Antarctic ozone changes. 


The models consistently simulate an enhancement of the heterogeneous chemical reactions in the 
GeoMIP experiments, with a direct significant depletion of stratospheric NO x , as reported in the 
literature after major volcanic eruptions [e.g. Johnston et al., 1992; Koike et a/., 1994]. Figure 5 shows 
the NO x (NO + NO2) mixing ratio changes for G4 during the central decade 2040-2050 (GISS-E2-R is 
not pictured in Fig. 5 because it does not include heterogeneous chemistry on the particle surfaces). All 
models calculate a decrease in NO x between 100 hPa and 5 hPa altitude in the tropics (20S-20N) and 
mid-latitudes (30S-50S and 30N-50N), of the order of 0.5-1 ppbv in the 10-50 hPa layer. Above 10 hPa 
the models show a small increase of NO x (0.2 to 0.5 ppbv at 2 hPa). This is most likely due to the fact 
that both O3 and NCb photolysis decrease as a consequence of the aerosol scattering (see ahead Fig. 7), 
so that (in addition to its enhanced heterogeneous chemical loss) less NO is available as sink for NO y 
through UV photolysis. This effect is less evident in GEOSCCM, which does not take into account the 
aerosol impact on photolysis. 

4.1 Changes in the vertical profile of ozone 

A comparison of O3 production/loss terms per chemical family in G4 respect to the base case has been 
made with the ULAQ model output. In order to highlight the direct impact of precursors changes on the 
ozone chemistry, the net chemical production term (P-LfCfi]) has been compared between the two 
simulations by keeping [O3] fixed at the RCP4.5 value, which means that a comparison of O3 
production (P) and loss frequency (L) terms is actually made (AP and [O3] AL are plotted in Fig. 6). 
The ULAQ-CCM calculates an increase in net ozone production rate in the tropics between 30 hPa and 
7 hPa for the decade 2040-2050 (on average +0.4-10 cm' s' ), and a decrease at all other tropical 
stratospheric altitudes (Fig. 6a). The significant NO x decrease discussed in Fig. 5 drives most of the 
chemical changes, by decreasing the O3 destruction from NO2 + O and the O3 production from NO + 
HO2 and NO + CH3O2. This results in a net increase of the ozone production term from the NO x cycles 
up to about 7 hPa (with a peak increase of + 1 .3 • 10 cm' s' at 20 hPa). 

Indirectly, the NO x change produced by the enhanced heterogeneous chemical reactions on sulfate 
aerosol SAD forces an increase of the O3 loss by FIO x , Cl x and Br x catalytic cycles (due to the 
decreasing loss of OH from OH + NO2 and the decreasing amounts of chlorine and bromine nitrates); a 
peak increase of -0.6- 10 cm' s' for these loss terms is calculated at 20 hPa. The perturbation to the 
O2 photolysis rate (see Fig. 7a) dominates below 50 hPa due to the UV screening effect by the aerosols; 
the ozone increase forced by the NO x cycle perturbation in the 7-30 hPa layer has a feedback on the O2 
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photolysis (Fig. 7a), thus decreasing the O3 production from JO2 in this same layer (-0.4- 10 cm' s' , 
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on average). Changes of ozone production/loss above 7 hPa (-0.2- 10 cm' s' , on average) result from 
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partially compensating effects of decreasing 0( P) from O3 photolysis (Fig. 7b) and increasing NO x 
(Fig. 5). 

Fig. 6b refers to the Antarctic springtime (September-October): in this case the largest perturbation 
takes place in Cl-Br catalytic cycles for ozone destruction (up to -1.0- 10 cm' s' at 50 hPa) and is due 
to the combination of increasing sulfate aerosol SAD, with decreasing NO x and increasing Cl x and Br x , 
diabatic cooling and enhanced formation of PSCs as a consequence of the local temperature decrease 
and the availability of additional sulfate aerosol condensation nuclei, with a positive feedback on ozone 
loss frequencies [ Tilmes et al., 2009]. 

Figure 8 shows the latitude-altitude distribution of annually averaged ozone net production changes. 
The layer of increased ozone net production between 7 and 30 hPa (up to +0.5-10 cm' s' ) is due to 
the suppression of the NO x cycle [Tie and Brasseur, 1995]. The layer of decreased ozone net 
production below 30 hPa (up to -0.5 - 1 0 cm' s' ) is due to a combination of reduced O2 photolysis and 
enhanced FIO x and Cl-Br cycles due to the NO x reduction. Negative changes above 7 hPa are produced 
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by coupled changes of NO x (Fig. 5) and of 0( P) from O3 photolysis (Fig. 7b). 

In the tropics, an important proportion of the ozone anomaly is caused by changes in upwelling due to 
perturbations of stratospheric heating rates [ Tilmes et al., 2009]. The heating from sulfate aerosols in 
the tropical middle stratosphere produces a temperature increase, which is mitigated by the ozone 
depletion in the same region. Anomalies of the global temperature at 70 hPa are shown in Fig. 9. GISS- 
E2-R simulates the smallest temperature anomaly. This is consistent with the smaller particle radius 
assumed in GISS-E2-R (see Table 2), which results in lower heating rates due to the smaller infrared 
absorption [Niemeier et al., 2011]. It is also consistent with larger ozone depletion, hence cooling, due 
to the larger change in UV scattering and O2 photolysis and to the absence of heterogeneous chemistry 
on sulfate aerosols . Figure 10 shows the simulated residual vertical velocity in the base case and its 
anomaly in G4 with respect to the base case. The models show a good consistency in the tropical 
upwelling and in the overall qualitative behavior of the G4 anomaly profile, with the magnitude of the 
largest relative increase of w* between 20 and 30 km altitude (50 and 10 hPa), ranging between 2 to 
5% for GISS-E2-R and MIROC-ESM-CHEM on one side, and 5 to 15% for ULAQ and GEOSCCM on 
the other side. 


Figure 1 1 shows the vertical profiles of ozone changes in G4 at the tropics and mid-latitudes, as 
simulated by ULAQ-CCM, MIROC-ESM-CHEM and GEOSCCM; GISS-E2-R results are discussed 
separately (see ahead Fig. 12), since this model does not include the effects of the enhanced sulfate 
aerosol SAD on the heterogeneous chemistry. In the tropical region (Fig. 11a and lib), ULAQ-CCM 
and MIROC-ESM-CFIEM predict a decrease of ozone between 100 hPa and 50 hPa mostly due to 
changes in CL photolysis (Fig. 7a). GEOSCCM simulates no change below 70 hPa, being the aerosol 
feedback on photolysis rates absent in this model. From 70 to 30 hPa the ozone negative changes due to 
increased upwelling (Fig. 10) contribute the total ozone anomaly together with the enhanced chemical 
loss in HO x , Cl x and Br x chemical cycles (Fig. 6a). The increased upwelling takes ozone poor air to 
ozone rich regions, creating the negative ozone anomaly that peaks at 30 hPa (-150 ppbv on average). 
Above 30 hPa the suppression of the NO x cycle dominates the chemical perturbation (about +200 ppbv 
at 10 hPa), with GEOSCCM showing a wider altitude range of the negative ozone anomaly, due to the 
larger increase of upwelling in this model (Fig. 10). At mid-latitudes (Fig. lie), where the effect on 
ozone of the increased upwelling disappears and the changes due to photolysis and heterogeneous 
chemistry dominate, GEOSCCM simulates smaller ozone anomalies with respect to ULAQ-CCM and 
MIROC-ESM-CCM. 

The ozone tropical changes in GISS-E2-R are compared to ULAQ-CCM (run in this case without 
heterogeneous chemistry on sulfate aerosols) in Fig. 12. The left panel shows the ozone anomalies in 
G3: the smaller aerosol particles in GISS-E2-R with respect to ULAQ-CCM (Table 2) scatter more UY 
radiation, increasing the O 2 photolysis (Fig. 7). This effect, however, is balanced by the lower AOD 
calculated by GISS-E2-R (Fig. 1), and results in an ozone anomaly vertical profile similar to ULAQ- 
CCM nhc. In G4 (Fig. 12, right panel), where the sulfate injection burden is fixed by the experiment 
design, the AOD simulated by GISS-E2-R and ULAQ-CCM is similar, and the net results is a five 
times larger ozone depletion in GISS-E2-R than in ULAQ-CCM. 

All models with heterogeneous chemistry simulate a significant ozone depletion in the Antarctic region 
(Fig. 13), due to a combination of increasing sulfate aerosol SAD (Fig. 13a) and enhanced fonnation of 
PSCs, produced in turn by local adiabatic and non-adiabatic cooling (Fig. 13b), the latter due to the 
feedback of photochemical ozone losses. ULAQ-CCM predicts the largest ozone depletion (Fig. 13c,d), 
both on annual basis and in springtime (-150 and -350 ppbv at 50 hPa, respectively). GEOSCCM 
simulates the smallest ozone depletion at 50 hPa during September-October, consistently with the 


lower temperature anomaly and the absence of feedback between PSCs occurrence and changes in 
sulfate aerosol concentration in the GEOSCCM. 

4.2 Column ozone changes and ozone radiative forcing 

Globally, the injection of geoengineering aerosol leads in all models to a reduction of the ozone column 
for the central decade of experiments G3 and G4 (2040-2050) (see Table 3). The stratospheric aerosol 
suppresses the ozone depleting NO x cycle, increasing ozone in the middle stratosphere, and enhances 
the C10 x and OH x cycles, decreasing ozone in the lower stratosphere [Tie and Brasseur, 1995]. 

Despite the constant stratospheric aerosol loading in G4, the magnitude of the geoengineering aerosol 
induced ozone depletion decreases in time (Fig. 14a) due to the decreasing atmospheric chlorine 
concentrations. ULAQ-CCM and MIROC-ESM-CHEM even simulates an increase in ozone starting 
from about 2050, when the ozone increase due to the suppression of the NO x cycle is no longer 
balanced by the decrease in total column due to C10 x and FIO x . 

On the other hand, the models simulates an increasing magnitude of ozone depletion in G3 (Fig 14b). If 
no heterogeneous chemistry is included in the simulations (ULAQ-CCM nhc and GISS-E2-R nhc), the 
increasing stratospheric aerosol burden leads to a lower ozone net production during the whole 
simulated period because of photolysis and temperature changes. If the heterogeneous chemistry is 
included (ULAQ-CCM, dashed line in Fig 14b), the depletion of the ozone column is reduced after 
2050, because of the increasing importance of NO x relative to C10 x . The net-adjusted tropopause 
radiative forcing (Tab. 3, Fig. 14d-e) is calculated off-line with the ULAQ-CCM radiative transfer code 
and its time behavior is closely correlated to that of the ozone column changes; it ranges between 0 and 
-0.1 W/m 2 (Fig. 14 d,e) except for the GISS-E2-R G4 case (not shown), where the O3 column 
perturbation is larger (see Table 3 and discussion in section 4.2) and consequently the tropopause net 
adjusted RF raises to about -0.23 W/m . These ozone RFs represents a rather small correction of the 
dominant direct forcing of geoengineering aerosols (about 2 to 10% of it). Table 3 suggests (as 
expected) that the largest contribution to the O3 RF comes the stratospheric temperature adjustment in 
the longwave range [see IPCC/TEAC, 2005]. 

Through the depletion of the ozone column, the GeoMIP aerosols cause an increase in the global UV-B 
radiation reaching the surface (Fig. 15a). Such increase is, however, mostly over-compensated (in the 
tropics) by the scattering of UV-B radiation by the aerosols (Fig. 15b). The net effect on the surface 
UV-radiation is not meridionally unifonn: while the tropics experience a decrease in UV-B radiation 


reaching the surface, the UV-B over Arctic and Antarctic polar regions increases by about 5% with 
respect to the base case in 2040-2050, following the ozone column perturbation. These findings are in 
agreement with Tilmes et al. [2012], who found a decrease of 5% erythemal UV in mid- and high 
latitudes for 2040 if the stratospheric halogen content from very short-lived halogen sources was 
considered. In these regions, the UV-B radiation changes are largely produced by the indirect ozone 
anomalies caused by geoengineering, but to a lesser extent also as direct consequence of the GeoMIP 
aerosols, which tend to increase the downward diffused radiation at large solar zenith angles [ Tsitas 
and Yung , 1996]. A well tested radiative transfer model (TUV) [ Madronich and Floke, 1998] has been 
used off-line for these UV-B calculations. 


5. Conclusions 

We have described the changes in ozone in simulations of GeoMIP experiments G3 and G4 in an 
ensemble of models. The results presented here are consistent with current understanding of the effects 
of different processes relevant to stratospheric ozone radiative, chemical, and dynamical effects. 

In the tropics, ozone changes are a complex combination of reductions due to changes in photolysis 
rates, increases due to suppression of the NO x cycle, reductions due to stratospheric heating and hence 
enhanced upwelling of ozone-poor air, and modulated reductions due to the addition of surfaces which 
serve as sites for heterogeneous chemistry. The net effect is a reduction of total column ozone in the 
tropics, which is enhanced in GISS-E2-R due to the lack of representation of heterogeneous chemistry, 
as well as a smaller aerosol size than was used in ULAQ-CCM. 

In the mid-latitudes, upwelling of ozone-poor air does not play a dominant role, so photolysis changes 
and heterogeneous chemistry dominate. GEOSCCM does not include photolysis changes, and as such, 
shows smaller changes in lower stratospheric ozone and total column. 

At the poles, particularly the Antarctic, cooling induced by photochemical ozone losses and circulation 
changes resulting from stratospheric heating enhances polar stratospheric cloud formation, thus 
promoting ozone depletion. The net effect is a reduction in polar ozone by over 5% in the majority of 
model simulations. 

The presence of a heterogeneous chemistry parameterization increases the importance of NO x relative 
to C10 x . As C10 x availability is reduced, depletion of the ozone column also lessens. In two models, 


the globally averaged ozone column increases after the year 2050. hi models which include 
heterogeneous chemistry, total ozone column changes are small (1-2 DU) and ozone radiative forcing 
is less than 0.1 W m' . As such, with the exception of polar regions, which show significant increases 
in UV-B, models predict small impacts on total column stratospheric ozone as simulated in the two 
experiments presented here. 

These results are potentially heavily dependent upon the experimental design. ULAQ-CCM is the only 
model included in this study which can represent aerosol microphysical growth; it shows a substantial 
increase in sulfate aerosol effective radius, consistent with past studies [Heckendorn et al., 2009; 
English et al., 2012]. A larger aerosol size reduces the radiative efficiency of the aerosols, requiring 
much more sulfate aerosol mass to achieve the desired radiative forcing [ Pierce et al., 2010]. 
Therefore, the effects on ozone described here may be magnified, depending upon the required mass of 
SO 2 to meet a climate goal. Moreover, these results are likely specific to continuous tropical injections. 
Modulating the latitudinal and temporal distribution of geo engineering could have different climate 
effects [MacMartin et al., 2013]. 

The results presented here would benefit from further uniformization of the experiment design. 
Although we have identified particular processes that impact the results, we have not been able to 
isolate the effects of the processes. For example, all models used different aerosol size distributions. It 
is known that larger particles have a greater infrared effect but a smaller surface area, so 
inhomogeneous representations of the aerosol size distribution confound changes in ozone due to the 
radiative or dynamical effects and changes due to heterogeneous chemistry. Obtaining robust features 
in a more controlled setting could be a useful follow-up study. 

The results presented here show only one aspect of geoengineering with stratospheric aerosols when 
conducted in a certain way. We do not advocate deployment of geoengineering, nor do we suggest 
ways in which it would be performed. Such decisions should be reserved for legitimate governance 
structures. Moreover, decisions regarding geo engineering should not be made based solely on physical 
science studies. Social, political, economic, and ethical perspectives also have an important role in 
informing decisions about geoengineering. 
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Tables 


Table 1. Summary of model runs (nhc = no heterogeneous chemistry on sulfate aerosols). 



RCP4.5 

G3 

G3-nhc 

G4 

G4-nhc 

ULAQ-CCM 

2 

2 

2 

2 

2 

GISS-E2-R 

3 


3 


3 

M1ROC-ESM-CHEM 

1 



4 


GEOSCCM 

2 



2 



Table 2. Stratospheric aerosol effective radius (pm) at 20 km over the tropics (2040-2050). For the 
observations, an average over 1999-2000 is used for the unperturbed background and over April 1992- 
March 1993 for a volcanic perturbation (i.e. Pinatubo) comparable in magnitude to G4 (in terms of 
average stratospheric mass burden of sulfate; see text). Background aerosol is not included in RCP4.5 
simulations with GEOSCCM (see Section 2.4 for details) 



RCP4.5 

G3 

G4 

ULAQ-CCM 

0.19 

0.48 

0.61 

GISS-E2-R 

0.15 

0.35 

0.35 

M1ROC-ESM-CHEM 

0.24 


0.24 

GEOSCCM 

* 


0.61 

SAGE-II 

0 . 22 ± 0.02 


0.59±0.02 


Table 3. O 3 column changes with respect to base case (DU); shortwave, longwave, adjusted longwave 
and net adjusted tropopause radiative forcing (mW/mj of ozone; surface UV-B change (percent) due 
to ozone changes. Annually and globally averaged values for the whole decade 2040-2050. 



0 3 col 
(DU) 

RF-O 3 SW 
(mW/m 2 ) 

RF-O 3 LW 
(mW/m 2 ) 

RF-O 3 LWadj 
(mW/m 2 ) 

RF-O 3 NET 

(mW/m 2 ) 

Surface UV-B 
change (%) 

G4 - RCP4.5 
ULAQ-CCM 

- 1.1 

2.0 

-16.5 

-34.1 

-32.1 

0.06 

G4 - RCP4.5 
MIROC-ESM-CHEM 

- 1.1 

8.8 

- 11.1 

-44.4 

-35.6 

0.22 

G4 - RCP4.5 
GEOSCCM 

- 2.1 

9.8 

-19.1 

-37.4 

-27.6 

0.51 

G3 - RCP4.5 
ULAQ-CCM 

- 2.8 

4.4 

-58.9 

-96.8 

-92.4 

1.16 

G4 nhc - RCP4.5 
ULAQ-CCM 

-1.9 

21.6 

-32.8 

-64.5 

-42.9 

0.88 

G4 nhc - RCP4.5 
GISS-E2-R 

-9.7 

127 

-23.6 

-361 

-234 

3.23 

G3 nhc - RCP4.5 
ULAQ-CCM 

-1.9 

22.7 

-29.5 

-60.9 

-38.2 

0.81 


G3 nhc - RCP4.5 
GISS-E2-R 


31.7 


-7.5 


-61.8 


-30.1 


0.69 


Figure captions 

Fig. labcde. Time series of globally averaged stratospheric AOD at 0.55 pm in the models: 1960-2000 
(top panel (a)) and G3, G4 perturbed cases in 2020-2070 (mid panels (b) and (c) respectively). Bottom 
panels (d), (e) show the zonally averaged AOD change with respect to RCP4.5, in G3 and G4, 
respectively. Line styles for the models are specified in the legend. 

Fig. 2ab. Vertical profiles of aerosol extinction ratio between 0.525 pm and 1.020 pm wavelength 
channels in the tropics (20S-20N), as calculated in the ULAQ-CCM for case G4 (solid line, 2040-2050 
average), for the Pinatubo case (dashed line, April 1992 - March 1993), and for background conditions 
(dash-dotted line, 1996-2000), compared to SAGE II retrievals (triangles and asterisks for 1992/1993 
and 1999/2000 conditions, respectively) (panel a). Panel (b): Aerosol size distribution (dn/dlogr, cm' ) 
calculated in the ULAQ model and averaged over the tropics and in the 30-70 hPa vertical layer (2040- 
2050). Solid line is for RCP4.5; dashed line is for case G4. 


2 3 

Fig. 3abc. Sulfate aerosol surface area density (jam /cm ) calculated in the models (ULAQ-CCM, 
MIROC-ESM-CHEM, GEOSCCM) for background conditions (1996-2000 average) and compared to 
SAGE-II data (asterisks). Panels (a,b,c) refer to tropics (20S-20N), mid-latitudes (30S-50S, 30N-50N) 
and polar regions (60S-90S, 60N-90N), respectively. 

Fig. 4abc. As in Fig. 3, but for but for case G4 (2040-2050 average). An indirect comparison is made 
with SAGE-II data for post-Pinatubo conditions (April 1992 - March 1993 average; triangles). The 
average sulfate mass loading above the tropopause is comparable in the latter case with respect to G4. 
The stratospheric mass burden comparison in these two cases is made using averages from ULAQ- 
CCM and GEOSCCM time-dependent simulations (4.83 ± 0.06 Tg-S for April 1992 - March 1993 
post-Pinatubo conditions and 4.7 ± 0.6 Tg-S for G4). 

Fig. 5ab. NO+NO 2 mixing ratio changes, as calculated in ULAQ-CCM, MIROC-ESM-CHEM and 
GEOSCCM (2040-2050). Panel (a) is for the tropics; panel (b) for mid-latitudes. 

Fig. 6ab. ULAQ-CCM calculated changes of net ozone production terms (10 mol cm' s' ) in G4 with 
respect to RCP4.5 (2040-2050); see legend for the different curves. Panel (a) is for the annual mean 
over the tropics (20S-20N); panel (b) is for the October mean over south polar latitudes (65S-90S). 

Fig. 7abc. Calculated percent changes of O2 and O3 — > 0( 3 P) photodissociation coefficients in panels 
(a) and (b), respectively, for G4 with respect to base case (tropics, spring equinox). A fixed O3 vertical 
profile (RCP4.5) is used in all cases, in order to evaluate the photodissociation sensitivity to the aerosol 
perturbation. Only the solid-circle line in panel (a) shows the JO 2 perturbation (in the ULAQ-CCM) 
with both aerosol and O3 changes taken into account in G4. 

Fig. 8. As in Fig. 6b, but for zonally averaged values of the total net production. 

Fig. 9ab. Panel (a): time series of calculated global temperature changes at 50 hPa, for G4 with respect 
to base case from 2020 to 2070 (K); solid, dashed, dash-dotted and dotted-circles line are for ULAQ- 
CCM, GISS-E2-R MIROC-ESM-CHEM and GEOSCCM, respectively. Panel (b): as above, but for 
ULAQ-CCM temperature anomalies from 1960 to 2000 relative to the 1980-1989 average compared to 
observations. Solid/dashed lines are with/without volcanic aerosol stratospheric heating, respectively. 
Observations are from sondes and ERA-40 reanalysis (solid line with diamonds) (Eyring et ah, 2006). 

Fig. lOab. Residual vertical velocity w* in the tropical stratosphere (mm/s) (average 15S-15N, 2040- 
2050) for baseline conditions in panel (a) and percent changes G4-BC in panel (b). ULAQ-CCM, 


GISS-E2-R, MIROC-ESM-CHEM and GEOSCCM and results are shown with solid, dashed, dash- 
dotted and dotted-circles line, respectively. 

Fig. 1 labc. Calculated ozone profile changes in G4 with respect to base case (2040-2050). Panel (a) is 
for tropical mixing ratio (ppbv); panel (b) for tropical number concentration (10 cm' ); panel (c) for 
mid-latitude number concentration (10 ~ cm' ). 

Fig. 12ab. As in Fig. 12b, but for G3 and G4 in panels (a) and (b), respectively, both without 
heterogeneous chemistry on sulfuric acid aerosols. 

Fig. 13abcd. G4-RCP4.5 changes averaged over the South Polar region (65S-90S) for 2040-2050. Top 
panels (a,b) are for aerosol surface area density (pm /cm ) and temperature (K), respectively (annual 
mean). Bottom panels (c,d) are for ozone mixing ratio (ppbv), annual and September-October mean 
values, respectively. Solid, dash-dotted and dotted-circles line are for ULAQ-CCM, MIROC-ESM- 
CFIEM and GEOSCCM, respectively. 

Fig. 14abcde. Top panels (a, b): time series of global O3 column changes in G4, G3, respectively, with 
respect to the base case (DU), averaged per decade. Mid panel (c): zonally and time averaged column 
ozone changes (2040-2050) with respect to RCP4.5. Bottom panels (d, e): as in panels (a, b), but for 
net adjusted tropopause RF (mW/mj. Line styles for the models are specified in the legend. 

Fig. 15ab. Zonally averaged UV-B percent changes at the surface, due to ozone and aerosol 
perturbations in G3, G4, with respect to RCP4.5 (average 2040-2050). Results for the different models 
are shown in top panel (a), with line styles as specified in the legend. A model average of the results is 
shown in the bottom panel (b), with dashed line for UV-B changes due to ozone perturbations only, 
dash-dotted line for aerosol perturbations only, solid line for the total. 
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Fig. labcde. Time series of stratospheric AOD at 0.55 pm in the models: 1960-2000 (top panel (a)) 
and G3, G4 perturbed cases in 2020-2070 (mid panels (b) and (c) respectively). Bottom panels 
(d), (e) show the AOD change as a function of latitude, with respect to RCP4.5, in G3 and G4, 
respectively. Line styles for the models are specified in the legend. 
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Fig. 2ab. Vertical profiles of aerosol extinction ratio between 0.525 pm and 1.020 pm wavelength 
channels in the tropics (20S-20N), as calculated in the ULAQ-CCM for case G4 (solid line, 2040- 
2050 average), for the Pinatubo case (dashed line, May 1992 - April 1993), and for background 
conditions (dash-dotted line, 1999-2000), compared to SAGE II retrievals (triangles and asterisks 
for 1992/1993 and 1999/2000 conditions, respectively) (panel a). Panel (b): Aerosol size 
distribution (dn/dlogr, cm' 3 ) calculated in the ULAQ model and averaged over the tropics and in 
the 30-70 hPa vertical layer (2040-2050). Solid line is for RCP4.5; dashed line is for case G4. 
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Fig. 3abc. Sulfate aerosol surface area density (pm /cm ) calculated in the models (ULAQ-CCM, 
MIROC-ESM-CHEM, GEOSCCM) for background conditions (1996-2000 average) and 
compared to SAGE-II data (asterisks). Panels (a,b,c) refer to tropics (20S-20N), mid-latitudes 
(30S-50S, 30N-50N) and polar regions (60S-90S, 60N-90N), respectively. 
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Fig. 4abc. As in Fig. 3, but for case G4 (2040-2050 average). An indirect comparison is made with 
SAGE-II data for post-Pinatubo conditions (April 1992 - March 1993 average; triangles), when 
the average sulfate mass loading above the tropopause is comparable to G4 (4.83 ± 0.06 Tg-S for 
April 1992 - March 1993 post-Pinatubo conditions and 4.7 ± 0.6 Tg-S for G4, average of ULAQ- 
CCM and GEOSCCM time-dipendent simulations). 



Fig. 5ab. NO+NO 2 mixing ratio changes, as calculated in ULAQ-CCM, MIROC-ESM-CHEM 
and GEOSCCM (2040-2050). Panel (a) is for the tropics (20S-20N); panel (b) for mid-latitudes 
(30S-50S, 30N-50N). 



Fig. 6ab. ULAQ-CCM calculated changes of net ozone production terms (10 s mol cm' 3 s' 1 ) in G4 
with respect to RCP4.5 (2040-2050); see legend for the different curves. Panel (a) is for the 
annual mean over the tropics (20S-20N); panel (b) is for the October mean over south polar 
latitudes (65S-90S). 



Fig. 7ab. Calculated percent changes of O2 (panel a) and O3 — > 0( 3 P) (panel b) photodissociation 
coefficients for G4 with respect to base case (tropics, spring equinox). A fixed O 3 vertical profile 
(RCP4.5) is used in all cases except for the solid-circle line, in order to evaluate the 
photodissociation sensitivity to the aerosol perturbation. The solid-circle line in panel (a) shows 
the J0 2 perturbation (in the ULAQ-CCM) with both aerosol and O 3 changes taken into account 
in G4. 
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Fig. 8. As in Fig. 6b, but for zonally averaged values of the total net production. 
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Fig. 9ab. Panel (a): time series of calculated global temperature changes at 50 hPa, for G4 with 
respect to base case from 2020 to 2070 (K); solid, dashed, dash-dotted and dotted-circles line are 
for ULAQ-CCM, GISS-E2-R MIROC-ESM-CHEM and GEOSCCM, respectively. Panel (b): as 
above, but for ULAQ-CCM temperature anomalies from 1960 to 2000 relative to the 1980-1989 
average compared to observations. Solid/dashed lines are with/without stratospheric heating 
from volcanic aerosol, respectively. Observations are from sondes and ERA-40 reanalysis (solid 
line with diamonds) (Eyring et al., 2006). 
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Fig. lOab. Residual vertical velocity w* in the tropical stratosphere (mm/s) (average 15S-15N, 
2040-2050) for baseline conditions in panel (a) and percent changes G4-BC in panel (b). ULAQ- 
CCM, GISS-E2-R, MIROC-ESM-CHEM and GEOSCCM and results are shown with solid, 
dashed, dash-dotted and dotted-circles line, respectively. 
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Fig. llabc. Calculated ozone profile changes in G4 with respect to base case (2040-2050). Panel 

12 3 

(a) is for tropical mixing ratio (ppbv); panel (b) for tropical number concentration (10 cm' ); 

12 3 

panel (c) for mid-latitude number concentration (10 cm' ). 


Ozone profile changes (2040-2050) 




j a1 2 -3 j a 12 -3 

10 cm 10 cm 


Fig. 12ab. As in Fig. 12b, but without heterogeneous chemistry on sulfuric acid aerosols for G3 
and G4 in panels (a) and (b), respectively. 




Fig. 13abcd. G4-RCP4.5 changes averaged over the South Polar region (65S-90S) for 2040-2050. 
Top panels (a,b) are for aerosol surface area density (pm /cm ) and temperature (K), respectively 
(annual mean). Bottom panels (c,d) are for ozone mixing ratio (ppbv), annual and September- 
October mean values, respectively. Solid, dash-dotted and dotted-circles line are for ULAQ- 
CCM, MIROC-ESM-CHEM and GEOSCCM, respectively. 
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Fig. 14abcde. Top panels (a, b): time series of global 0 3 column changes in G4, G3, respectively, 
with respect to the base case (DU), averaged per decade. Mid panel (c): zonally and time 
averaged column ozone changes (2040-2050) with respect to RCP4.5. Bottom panels (d, e): as in 
panels (a, b), but for net adjusted tropopause RF (mW/m 2 ). Line styles for the models are 
specified in the legend. 
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Fig. 15ab. Zonally averaged UV-B percent changes at the surface, due to ozone and aerosol 
perturbations in G3, G4, with respect to RCP4.5 (average 2040-2050). Results for the different 
models are shown in top panel (a), with line styles as specified in the legend. A model average of 
the results is shown in the bottom panel (b), with dashed line for UV-B changes due to ozone 
perturbations only, dash-dotted line for aerosol perturbations only, solid line for the total. 


